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ABSTRACT 

Wc investigate the relations among the stellar continuum-subtracted 8 /im polycyclic 
aromatic hydrocarbon (PAH 8 /im) emission, 24 /im hot dust emission, and 160 //m 
cold dust emission in fifteen nearby face-on spiral galaxies in the Spitzer Infrared 
Nearby Galaxies Survey sample. The relation between PAH 8 and 24 /im emission 
measured in ~ 2 kpc regions is found to exhibit a significant amount of scatter, and 
strong spatial variations are observed in the (PAH 8 /im)/24 /im surface brightness 
ratio. In particular, the (PAH 8 /im)/24 /im surface brightness ratio is observed to 
be high in the diffuse interstellar medium and low in bright star-forming regions and 
other locations with high 24 fim surface brightness. PAH 8 /im emission is found to 
be well-correlated with 160 /im emission on spatial scales of ~ 2 kpc, and the (PAH 
8 /im)/160 /im surface brightness ratio is generally observed to increase as the 160 /im 
surface brightness increases. These results suggest that the PAHs are associated with 
the diffuse, cold dust that produces most of the 160 /im emission in these galaxies, 
and the variations in the (PAH 8 /im)/160 /im ratio may generally be indicative of 
either the intensity or the spectrum of the interstellar radiation field that is heating 
both the PAHs and the diffuse interstellar dust. 
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1 INTRODUCTION 

Observations of nearby spiral galaxies with the Infrared 
Space Observatory (ISO) demonstrated that polycyclic aro- 
matic hydrocarbon (PAH) spectral feature emission at 
~ 8 /xm is closely associated with hot dust emission at 
15 /im, and the 8 /im emission was also sho wn to be cor- 
related with other star formation tracers (|Roussel et al.l 
l200ll : iForster Schreiber et al. I l2004 l h_ However, observ at ions 
with the Spitzer Space Telescope dWerner et al.ll2004l ) have 
demonstrated that the correlation of PAH emission at 
8 /zm to hot dust emission at 24 /im was not neces- 



sarily a one - to-on e correlation within indivi dual galaxies. 
iHelou et all i|2004l h iBendo et all i|2006] h and iGordon et ail 
l|2008h demonstrated that, on scales of hundreds of parsecs, 
continuum-subtracted PAH emission at 8 /tm (henceforth re- 
ferred to as PAH 8 )im emission) does not accurately trace 
24 fim emission. In particular, PAH 8 /im emission seems 
to appear in shell-like features around star-forming regions, 
while 24 /im emissi on peaks within star- forming regions. Us- 
ing images of M51, ICalzetti et al] (|2005) demonstrated that 
PAH 8 /im emission compared to other tracers of star for- 
mation, including 24 fim emission, is disproportionately low 
within star-forming regions and disproportionately high in 
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diffuse regions. The results from lCalzetti et al. (2007) show 
that the PAH 8 /xm/24 /mi emission ratio measured in HII 
regions is partly dependent on the aperture used for measur- 
ing the flux densities, which also suggest that a significant 
fraction of PAH 8 /xm emission originates from outside the 
HII regions. iThilker et all <|2007l ) also demonstrated that the 
PAH 8 /im/24 /xm ratio of HII regions in NGC 7331 may de- 
pend on infrared, Ha, or ultraviolet surface brightness after 
corrections for dust extinction have been applied, which im- 
plied that PAH 8 /xm emission is suppr e ssed w i thin strong 
star- fo r ming regions. lEngelbracht et ail (I2005T) . iDale et al.l 
l|2005l ). IDraine et all (|2007h . and lEngelbracht et all (|2008T l 
have demonstrated that the global 8 /xm/24 /xm emission 
ratio may decrease in regions with lower metallicity, and 
ICalzetti et al.l (|2007h also demonstrated that this ratio for 
HII regions in other galaxies depends on metallicity. Al- 
though these studies based on Spitzer data have shown ex- 
plicitly or implicitly that the relation between 8 /im PAH 
emission and 24 /xm emission on <jl kpc scales is not a one- 
to-one relation, no study has yet explored variation in the 
(PAH 8 /xm)/24 /xm ratio between diffuse and star-forming 
regions within a broad range of spiral galaxies. 

While most research has focused on comparisons be- 
tween PAH emission and hot dust emission, some additional 
studies have compared PAH emission to cold dust emis- 
sion at wavel e ngths longe r than 100 / xm. U sing ISO data, 
iMattila et ail (119991 ^ and lHaas et al.l l|2002l ) demonstrated 
that PAH em ission at 8 / xm is correlated with large grain 
emission, and lHaas et alj (|2002T ) even argued that the cor- 
relation between 8 /xm emission and cold dust emission at 
850 (im was much stronger than the co rrelation between 8 
and 15 /xm emission. Using Spitzer data. lBendo et all (120061 ) 
demonstrated that PAH 8 /xm emission was correlated with 
160 /xm cold dust emission in NGC 4631, but the ratio of 
(PAH 8 /xm)/160 /xm emission varied. Although the (PAH 
8 /xm)/160 /xm ratio app eared related t o the infrared surface 
brightness in NGC 4631, lBendo et all ([20061 ) used bright in- 
frared sources in the outer disc to argue that the variations 
in the (PAH 8 /xm)/160 /xm ratio are instead dependent 
on radius. Given that the ratio of PAH emission to dust 
emission at longer wavelengths or the ra tio of PAHs to to- 
tal dust mass may vary with metallicity (lEngelbracht et all 
20051; IDale et al.ll2005l;rSmith et al.ll2007l ; IDraine et al.ll2007l ; 
Engelbracht et al.ll20081 ) and that metallicity is gener ally ex 
pected to vary with radius w i thin spiral galaxies (e.g. SearL 
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Zaritskv et~ail Il994l ; Ivan Zee et"ail Il998l ; iPilvugin et~ail 
20041 . J. Moustakas et al. 2008, in preparation), the re- 
sults from NGC 4631 suggested that the variations in (PAH 
8 /xm)/160 /xm ratio within that galaxy could be related to 
metallicity. Unfortunately, NGC 4631 is viewed edge-on, so 
local variations in the observed (PAH 8 /xm)/160 /xm ra- 
tio that are not dependent on radius may be suppressed by 
the line-of-sight integrations through the disc. Variations in 
the (PAH 8 /xm)/160 /xm ratio and its dependence on both 
infrared surface brightness and radius need to be studied 
further using spiral galaxies with orientations closer to face- 
on. 

In this paper, we examine the relation of PAH 8 /xm 
emission to both 24 /xm hot dust emission and 160 /xm cold 
dust emission in a set of fifteen face-on spiral galaxies ob- 
served by Spitzer as part of the Spitzer Infrared Nearby 



Gala xies Survey (SINGS) legacy project (|Kennicutt et al] 
2003). The basic goal is to understand how PAH 8 /xm emis- 
sion is related to dust emission at other wavelengths. Sec- 
tion [5] provides basic information on the wave bands used 
in the analysis, the observations and data reductions, the 
sample, and the preparation of the data for the analysis. 
The comparison between PAH 8 and 24 /xm emission is pre- 
sented in Section |31 and the comparison between PAH 8 and 
160 /xm emission is presented in Section [4] The analysis is 
then followed by a discussion in Section [S] and conclusions 
in Section [S] 



2 DATA 

2.1 Wave band information 

Here we present background information on the 8 /xm, 
24 /xm, and 160 /xm bands used for the analysis in this paper. 
This information also includes some caveats that should be 
considered when interpreting the results from these data. 

Channel 4 of the Infrared Array Camera (IRAC; 
iFazio et al.ll2004T l covers a region centered on 8.0 /xm that 
includes the 7.7 /xm PAH feature as well as some stellar con - 
tinuum and some hot dust emission fe.g. lSmith et al.ll2007n . 
For the analysis here, the stellar continuum is subtracted 
using 3.6 /xm data from channel 1 of IRAC. Although the 
PAH emission features are very prominent in this 8 /xm band, 
thermal emission from larger grains may contribute signifi- 
cantly to emissi on in this band in regions with very intense 
radiation fields (|Draine fc Lill2007h . 

The 24 /xm detector of the Multiband Imaging Pho- 
tometer for Spitzer (MIPS; lRieke et al.lliooi ') mostly detects 
hot ( J>100 K) dust emission within nearby galaxies. In envi- 
ronments with low radiation fields, this dust emission may 
originate from mostly transiently-heated small grains, but 
in regions with high radiations fields such as star-forming 
regions, the emission consists of mostly thermal emission 
from grains at equilibrium temperat ures of JjlOO K (e.g. 
iLi fc Drainel 1200 ll ; IDraine fc Lil 12007'). The relative contri- 
bution of dust at JjlOO K is expected to be strongly de- 
pendent on the strength of the illuminating radiation field. 
Consequently, the 24 /xm band is predicted to increase more 
rapidly than other infrared bands as the illuminating radi- 
ation field increases, as has been shown with ma ny physical 
and semi-empirical models of dust emission (e.g. IDale et all 
l200ll : ILi fc Drainel l200ll ; IDraine fc Lill2007T l. Multiple stud- 
ies with Spitzer data have shown that the 24 /xm emission 
from point-like sources corresponds to HII regions visible in 
optical and ultraviolet wave bands, so the 24 /xm band may 
be used to measure star formation activity (|Prescott et al.l 
120071 ; ICalzetti et al.ll2007h . Extended 24 /xm emission, how- 
ever, may originate from outside of star-forming regions in 
the diffuse interstellar medium. 

The 160 /xm MIPS detector mainly traces cold (~ 20 K) 
dust emission that may be asso ciated with cirrus d ust in the 
diffuse interstellar medium (e.g. IDraine et aill2007f ). Because 
the 160 /xm band samples dust emission near the peak of 
the spectral energy d istribution in the galaxies in this sam- 
ple l|Dale et al.ll2005l) . the 160 /xm band should be tightly 
correlated with the total infrared (TIR) luminosity. This is 
explained in the following proof. For thermal emission mod- 
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ified by an emissivity function that scales as v@ , the total 
energy emitted in a single wave band iyLJ) should scale as 



vL v oc 



where 



,(4+/3) 



X = 



h v 
kf' 



(1) 



(2) 



The total integrated thermal emission Ltotal will scale ac- 
cording to 

iwccT< 4+ ». (3) 

The ratio of vL y to L to tai will be approximately constant 
(e.g. independent of grain temperature) when it is close to 
a maximum, which occurs when 



e x — 1 

Since e x >> 1, then 
„ 1 hv 



(4) 



(5) 



4 + /3 k 

Assuming that f3 equals 2 (|Li fc Dr ainc 2001), this condition 
is met at 160 /im when the dust temperature is ~ 15 K. Be- 
cause the peak of the spectral energy distribution is actually 
at a shorter wavelength than 160 /jm, the ratio of 160 /im 
emission to total infrared emission should decrease slightly 
as the illuminating radiation field increases. 

2.2 Observations and data reduction 

We use the 3.6, 8.0, 24, 70, and 160 /im data taken with 
Spitzer as part of SINGS. The 3.6 and 8.0 /im observations 
were performed with IRAC. The observations for each ob- 
ject consisted of a series of 5 arcmin x 5 arcmin individ- 
ual frames taken in either a mosaic or a single field dither 
pattern. The 24, 70, and 160 /jm observations performed 
with MIPS are composed of two scan maps for each tar- 
get. Each object was observed twice in each wave band to 
identify and remove transient phenomena, particularly as- 
teroids. The full-width half-maxima (FWHM) of the point 
spread f unctions (PSFs), as stated in the Spitzer Observer's 
Manual (|Spitzer Science Centerll2006! fl. are 1.7, 2.0, 6, and 
38 arcsec at 3.6, 8.0, 24, and 160 /im, respectively. Details 
on the observations can be foun d in the documenta tion for 
the SINGS fourth data delivery (|SINGS Teaml l2006fl. 

The IRAC data were processed using the SINGS IRAC 
pipeline, which combines multiple frames of data using a 
dr izzle technique. A d escription of the technique is presented 
in iRegan et alj l|2006l ). The final images may contain resid- 
ual background emission from the telescope or sky that is 
subtracted during the analysis. The MIPS data were pro- 
cessed using the M IPS Data Analysis Tools version 3.06 
l|Gordon et all [2005). Additional software was used to re- 
move zodiacal light emission and improve the flatfielding in 
the 24 /im data and to remove short-term variations in the 
background signal (commonly referred to as drift) in the 70 
and 160 /im data. Any additional background offset left in 

6 http:/ /ssc. spitzer. caltech.edu/documents/som/ 

7 Available at http:/ /ssc. spitzer. caltech.edu/legacy/singshistory.html 




the final images was measured in regions outside the opti- 
cal discs and subtr acted. Additional details are presented in 
iBendo et al. (2006). Full details on the data processing are 
also available in the SINGS docu mentation for the fourth 
data delivery (|SINGS Teamll2006h . 

2.3 Sample selection 

To perform this analysis, we need to resolve substructures 
in the 160 /im images, which have PSFs with FWHM of 
38 arcsec. We therefore limit the sample to spiral galaxies in 
SINGS where the major axe s of the D25 isophote specified by 
Ide Vaucouleurs et al.l (|l99ll ) are larger than 5 arcmin. Since 
we want to be able to distinguish between radial colour vari- 
ations and colour variations related to the presence of sub- 
structures and since such substructures are difficult to study 
in edge-on galaxies, we only use galaxies that are inclined 
less than ~ 60deg. The inclinations are calculated using 



(6) 



The value q is the observed (projected) minor-to-major axis 
ratio. The value q is the intrinsic optical axial ratio (the ra- 
tio of the unprojected optical axis perpendicular the plane 
of the galaxy to the diameter of t he disc), wh ich is equiva- 
lent to 0.20 for most disc galaxies l|Tullvll 19981 ). Because the 
optical disc of NGC 5194 may be distorted by its interac- 
tion with NGC 5195, its inclination is not calculated using 
this equation. I nstead, the inclination as w ell as the position 
angle given bv lGarci'a-Gomez et al.l <|2002h are used for the 
analysis. 

Six of the SINGS galaxies that meet the above criteria 
are unsuitable for the analysis. NGC 1512 and NGC 4826 
are not used because only the central regions were detected 
at the 5<t levels in all of the convolved maps (described in the 
next section). The 8.0 /im images of NGC 1097, NGC 1566, 
and NGC 4736 are heavily affected by muxbleed artefacts 
(artificially bright columns of pixels associated with high- 
surface brightness sources) that cross over significant frac- 
tions of the optical discs, so those data are not usable for 
this analysis. Two very bright foreground stars in the 3.6 and 
8.0 /xm images of NGC 3621 cause problems in the analysis, 
so NGC 3621 needs to be excluded from the sample as well. 
The other 15 galaxies that meet the above criteria, which 
are roughly uniformly distributed between Hubble types Sab 
and Sd, are listed in Table [T] along with information on the 
galaxies' morphologies, optical axes, distances, nuclear spec- 
tral types, and nebular oxygen abundances (12+log(0/H), 
which is treated as representative of the global metallicities 
of the galaxies). 



2.4 Data preparation 

Many Spitzer studies of infrared colour variations within 
individual galaxies have relied on flux densities measured 
within discrete subregions that are chosen either by eye 
or by source identification software. However, this selection 
process may be biased. Some subregions within the galaxy 
might be excluded from the analysis, especially if the re- 
gions are selected by eye, and region selection may also be 
biased towards regions that appear bright in a specific wave 
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Table 1. Basic Properties of the Sample Galaxies 



Name 


niiDDie 


Size of Optical 


Inclination^ 


Position 


Distance 


Distance 


Nuclear 


1 0M]f\cr(d /tT\ 




lype 


Disc (arcmin) a 


Anglc c 


^IVlpC ) 




Reference e 


lype J 




NGC 628 


SA(s)c 


10.5 x 9.5 


26° 


25° 


7.3 ± 1.4 


i 


SF 


8.33 ± 0.02 


NGC 925 


SAB(s)d 


10.5 x 5.9 


58° 


102° 


9.12 ±0.17 


2 


SF 


8.24 ±0.01 


NGC 2403 


SAB(s)cd 


21.9 x 12.3 


58° 


127° 


3.13 ±0.14 


2 


SF 


8.31 ±0.01 


NGC 3031 


SA(s)ab 


26.9 x 14.1 


60° 


157° 


3.55 ±0.15 


2 


AGN 


8.40 ±0.01 


NGC 3184 


SAB(rs)cd 


7.4 x 6.9 


22° 


135° 


11.1 ± 1.9 


3 


SF 


8.48 ±0.02 


NGC 3351 


SB(r)b 


7.4 x 5.0 


49° 


13° 


9.3 ±0.4 


2 


SF 


8.61 ±0.01 


NGC 3938 


SA(s)c 


5.4 x 4.9 


25° 


42°* 


13.4 ±2.3 


4 


SF 


8.44 ±0.18 


NGC 4254 


SA(s)c 


5.4 x 4.7 


30° 


68°* 


16.5 ±0.6 


5 


SF 


8.46 ±0.02 


NGC 4321 


SAB(s)bc 


7.4 x 6.3 


32° 


30° 


14.3 ±0.5 


2 


AGN 


8.50 ±0.04 


NGC 4579 


SAB(rs)b 


5.9 x 4.7 


38° 


95° 


16.5 ±0.6 


5 


AGN 


8.55 ±0.18 


NGC 4725 


SAB(r)ab pec 


10.7 x 7.6 


46° 


35° 


11.9 ±0.3 


2 


AGN 


8.40 ±0.05 


NGC 5055 


SA(rs)bc 


12.6 x 7.2 


57° 


105° 


7.8 ± 2.3 


4 


AGN 


8.42 ±0.04 


NGC 5194 


SA(s)bc pec 


11.2 x 6.9 


20° 


170° 


8.4 ±0.6 


6 


AGN 


8.54 ±0.01 


NGC 6946 


SAB(rs)cd 


11.5 x 9.8 


32° 


80°* 


6.8 ± 1.7 


7 


SF 


8.40 ±0.04 


NGC 7793 


SA(s)d 


9.3 x 6.3 


49° 


98° 


3.9 ±0.4 


8 


SF 


8.27 ±0.03 



a These data are taken from lde~Vaucoulcu rs et al ] jl99ll V The optical disc is the size of the D25 isophote. 

b The inclinations are calculate d using the dimensio ns of the galaxies and Equation \E\ An exception is made for NGC 5194, 
where the inclination fromlGarci'a-Gomc z et al . (2002) was used. 

c The position angles are measured in degrees from north through ea st. These data are taken f rom Ide V aucoulcu rs et al.l 
l|l99lh when given. The position angle for NGC 5194 was taken from Garci'a-Go mez et al ] j2002h . Position angles marked 
with an asterisk are measured in the 3.6 fim data using the IDL fit_ellipse program by D. Fanning (accessible from 
http://www.dfanning.com/documents/programs.html). These position angles have uncertainties of ~ 5°. 

d These distances are part of a compilation of distance for the complete SINGS sample will be presented in J. Moustakas et al. 

(2008, in preparation). 

e Distance refer ences - CD ISharina et al] lll996l) ; (2) [Freedman et alJ fcOOll ); (3 ) iLeonard et al.l [|2002ri; (4 ) iMasterj d2005Tl ; (5) 
iMei et al] (|2007h : (6) iFeldmeier et al] jl997h : (7) iKarachentsev et al.l l|2000h : (81 iKarachentsev et al] l|2003h . 

f The nucleus types are based on the analysis of the optical spectrophotometric properties by J. Moustakas et al. (2008, in 
preparation). The spectral type indicates whether the center of each object is dominated by star formation (SF) or an active 
galactic nucleus (AGN). See J. Moustakas et al. (2008, in preparation) for more details. 

9 These are the charact eristic nebular oxygen abundances presented by J. Moustakas et al. (2008, in preparation) using the 
iPilvugin fe Thuanl {2005) strong-line calibration. They are representative of the mean, luminosity-weighted oxygen abundance of 
each galaxy. Note that the uncertainties only include statistical measurement errors, and do not include any systematic uncer- 
tainties in the adopted abundance scale (s ee J. Moustakas et al. 2008, i n preparation for more details). Although J. Moustakas 
et al. also prese nt abundances b ased on the lKobulnickv fc Kewlevl l|2004h strong-line calibration, we adopt the abundances bas ed 
on the Pilyugin & Thuan (2005) calibration because they are comparable to the abundances used by Engclbracht e t al.l (|2008| ). 



band. To avoid this bias, we divide each galaxy into 45 arc- 
sec square regions and extract surface brightnesses for all 
usable regions within the optical discs of the galaxies. This 
approach samples the whole of the galactic disc and avoids 
any subjective biases. We choose to use 45 arcsec regions 
because it is an integer multiple of the pixels used in the 
original images, and it is larger than the 38 arcsec FWHM 
of the 160 /Urn images. Note that 45 arcsec corresponds to 
physical scales of ~ 0.7-3.6 kpc for the galaxies in this sam- 
ple. 

First, the data are convolved with kernels that match 
the PSFs of the images in the 3.6, 8, 24, and 70 (im bands to 
the PSF of the 160 [im data, which has a FW HM of 38 arc- 
sec. T he convolution kernels were created by (I Gordon et al.l 
2008) using the ratio of the Fourier transforms of the input 
and output PSFs, with high-frequency noise suppressed in 
the input PSFfl By doing this, it is possible to directly com- 
pare surface brightnesses measured in the same apertures in 
different wave bands without applying any aperture correc- 
tions. Next, the coordinate systems across all the wave bands 



are matched to each other using point-like sources (stars, 
background galaxies, or infrared-bright regions within the 
galaxies) or the centres of the galaxies as guides. The back- 
ground in the IRAC data was then subtracted from the im- 
ages. Following this, the data are rebinned into 45 arcsec 
square regions such that the coordinates of the rebinned 
pixels match across all wave bands. Because the centres of 
the galaxies are mapped into the corners of four pixels in the 
MIPS images, the centres of the galaxies fall at the corners 
of four 45 arcsec bins. We then extract surface brightnesses 
from all 45 arcsec regions within the optical discs of the 
galaxies. We exclude from the analysis regions not detected 
at the 3a level in one or more wave bands, regions contam- 
inated by emission from bright foreground stars (identified 
as unresolved sources in the unconvolved IRAC data with 
3.6 /im/8 /im surface brightness ratios J>5), and regions af- 
fected by muxbleed at 8 /im. Additionally, we exclude the 
centre of NGC 5195 from the analysis of NGC 5194 be- 
cause the dwarf galaxy causes confusion when interpreting 
the data on the spiral galaxy. To correct for the diffusion 
of light through the IRAC detector substrate, the 3.6 and 
8.0 fim data are mul tiplied by the "infi nite" aperture cor- 
rections described bv lReach et al] (|2005l ) (which should not 



8 The kernels are available at 

http: / / dirty, as. arizona. edu / ~kgordon / mips / conv_psfs / conv_psfs .html 
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Table 2. Uncertainties in the Surface Brightness Ratios Related 
to Matching Image Coordinate Systems 



Surface 


Uncertainty 


Brightness 




Ratio 




7"„(PAH 8 fim)/ 1^(24 pm) 


0.005 


A, (PAH 8 imw)/Iu(1G0 fim) 


0.0015 


^(PAH 8 /im)/J(TIR) 


0.007 


Iu(24 j um)//„(160 pn) 


0.005 


y/„(160 /tm)//(TIR) 


0.02 



be confused with the types of aperture corrections used for 
measuring the flux densities of unresolved sources). The cor- 
rection factors are 0.944 and 0.737 at 3.6 and 8.0 jim, re- 
spectively. Finally, we subtract the stellar continuum from 
the 8 and 24 /im surface brightnesses (in MJy sr" 1 ) using 

I V (PAH 8/jrrc) = I v (8fj,m) - 0.232J„ (3.6/zm) (7) 

/„(24/xm) = I v (24fj,m) - 0.032I„ (3. 6/xm), (8) 

which were derived by iHelou et all (|2004j ). For the galax- 
ies in this analysis, this is a correction of ~ 1 %-4% to the 
24 /im surface brightnesses but a correction of ~ 5%-25% 
to the 8 fim surface brightnesses. The stellar continuum- 
subtracted 8 fim data are referred to as the PAH 8 /jm emis- 
sion throughout this paper. While these stellar continuum- 
subtracted spectra may contain some thermal dust emission, 
most of the emission should originate fr om the PAH emissio n 
features, as has been demonstrated by ISmith et al.l (|2007h . 
We also calculated TIR surface brightnesses using 

I{TIR) = 0.95vI„(PAH8 fim) + 1. 15^(24 fim) 
+vl v (70 fan) + W„(160 fim) *- ' 

based on equation 22 from iDraine fc Lil (|2007h . 

In the analysis of these 45 arcsec regions, we assume 
that the background noise, which is measured in off-target 
regions in the images, is the only source of uncertainties 
that will affect this analysis. Uncertainties in the calibra- 
tion will only scale the data, which does not affect either 
the slopes of lines fit to the data in log space or the scatter 
around the best fit lines. Uncertainties in the mean back- 
ground value subtracted from the data should be negligi- 
ble. The uncertainties for the ratios of surface brightnesses 
in two wave bands also includes a term to account for un- 
certainties in matching the coordinate systems of the two 
bands. This term is estimated from the standard deviations 
of the central four pixels in galaxies with symmetric, bright, 
point-like nuclei. In galaxies with such nuclei, the central 
four pixels should each sample approximately one-quarter 
of the peak of the central PSF. They should have the same 
surface brightness ratios if the image coordinate systems of 
different images are matched properly, but in practice, the 
ratios vary by small amounts. These terms, which are given 
in Table[2]for the ratios used in this paper, add uncertainties 
of approximately 5%-10%. 

We also use maps of the (PAH 8 fim)/ 24 fim and (PAH 
8 /im)/160 fim surface brightness ratios as additional aids 
in the interpretation of the data. To make these ratio maps, 
we again use the convolution kernels of K. D. Gordon to 
match the PSFs of the 3.6 and 8.0 /im images to the PSF 



in the longest wave band. Hence, the 3.6, 8.0, and 24 /im 
images used to produce the (PAH 8 /im) /24 fim maps have 
PSFs with FWHM of 6 arcsec, and the 3.6, 8.0, and 160 fim 
images used to produce the (PAH 8 /im)/160 /jm maps have 
FWHM of 38 arcsec. The coordinate matching, IRAC back- 
ground subtraction, IRAC aperture correction, and stellar 
continuum subtraction steps applied to the surface bright- 
nesses measured in the 45 arcsec regions are also applied 
to the images used to make these ratio maps. Only pixels 
where the 24 fim data were detected at the 5cr level and 
where the PAH 8 and 160 /im data were detected at the 
10cr level are shown in the maps. These thresholds were se- 
lected to filter out background noise and artefacts, which 
is why the threshold differs among the three wave bands. 
Moreover, pixels above these thresholds roughly correspond 
to the parts of the disc that were binned into the 45 arcsec 
regions described above. Additionally, foreground stars and 
any bright artefacts that remained in the data were masked 
out in the final ratio maps. Images of the PAH 8 fim emis- 
sion, the (PAH 8 fim)/ 24 /mi surface brightness ratio, and 
the (PAH 8 /mi)/ 160 fim surface brightness ratio are shown 
in Figure [T] 



3 COMPARISONS OF PAH 8 AND 24 fM EMISSION 

The panels in the middle columns of Figure [T] show how the 
(PAH 8 /im)/24 fim surface brightness ratio varies within 
the discs of the galaxies. The resolution of the maps is 6 arc- 
sec, which allows for seeing fine details within the images. 
If the PAH 8 /jm surface brightness varies linearly with the 
24 fim surface brightness, then this quantity should be con- 
stant across the discs of these galaxies. Instead, the (PAH 
8 pm)/24 /xm ratio is seen to vary significantly. In partic- 
ular, the 24 /mi emission peaks more strongl y in point-like 
region s, wh ich, according to the results from ICalzetti et al.l 
(|2007t) and iPrescott et all (120071 ). are HII regions. In con- 
trast, the PAH 8 /im emission is stronger relative to the 
24 /im emission in the diffuse interstellar regions. This is 
most apparent when comparing the arm and interarm re- 
gions of the grand-design spiral galaxies in this sample, no- 
tably NGC 628, NGC 3031, and NGC 5194. Some HII re- 
gions well outside the centres of the galaxies have very low 
(PAH 8 /im) /24 fim ratios even compared to the rest of the 
optical discs of the galaxies, as can be seen with the HII re- 
gion in the northeast part of the disc in NGC 3184 and the 
HII region in the east part of the disc in NGC 3938. 

Plots comparing colour to surface brightness are used 
here and in the next section to study the correlation between 
two wave bands. In a first-order approximation, the PAH 8, 
24, and 160 jim surface brightnesses are all correlated with 
one another for these data. Hence, colour variations may 
be difficult to see or measure in plots directly comparing 
surface brightnesses in two wave bands. This is because the 
surface brightnesses vary by factors of ~ 100 but the ratios 
of the surface brightnesses vary by less than a factor of 10. 
However, plots comparing surface brightness ratios to sur- 
face brightnesses can more clearly show deviations in the 
relation between two wave bands, including systematic vari- 
ations in colour related to surface brightness. If a one-to-one 
correspondence exists between the PAH and 24 /tm surface 
brightnesses, then the slopes of the best fit lines in these 
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Figure 1. Images of the PAH 8 fim emission at the resolution of the IRAC data (left column), the (PAH 8 fim)/24 /an surface brightness 
ratio at the resolution of the 24 /an data (centre column), and the (PAH 8 fim) /WO fim surface brightness ratio at the resolution of 
the 160 (im data (right column) for the 15 galaxies studied in this analysis. The three maps for each galaxy are scaled to the same size. 
North is up and east is left in each map. The stellar continuum has been subtracted from the 8 and 24 fim bands to produce these maps. 
The contour levels in the ratio maps are chosen to show structure without showing excessive scatter from noise. The colour bars on this 
page give the values of the contour levels. The red colours in the ratio maps correspond to relatively weak PAH 8 /im emission, and 
the blue colours correspond to relatively strong PAH 8 /an emission. White regions in the ratio maps correspond to regions with low 
signal-to-noise ratios, regions with foreground stars, or regions strongly affected by artefacts in the data. In the PAH 8 fim image of NGC 
3351, strips of data to the northeast and southwest of the centre have been strongly affected by muxbleed and have been masked out 
in the ratio maps. The circle in the top right corner of the (PAH 8 /^m)/160 (im ratio maps shows the 38 arcsec FWHM of the 160 fim 
beam. 



plots would be close to 0, and the scatter in these plots 
would be small. Such plots would be biased towards pro- 
ducing relations with slopes of -1 in log-log space if the two 
wave bands were randomly distributed. However, since the 
flux densities in all wave bands we studied here are approx- 



imately proportional to each other, such biases will not be 
present. To measure the physical (intrinsic) scatter around 
the best fit line, we will subtract the sum of the square of 
the measurement uncertainties from the sum of the square 
of the residuals from the best fit line. This is given by 
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Figure 1 - continued 
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Figure 1 - continued 
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Figure 2. Plots of the (PAH 8 fim)/ 24 /an surface brightness ratios versus the 24 fim surface brightnesses for the 45 arcsec square regions 
measured in these galaxies. Note that the data are extracted from images with the same resolution as the 160 fim data. This was done 
so that resolution effects will not be a factor when comparing the data in this figure to the data in Figure [5] The grey lines are the best 
fitting lines for the relations in each plot; slopes and intrinsic scatters for these fits are given in Table [3] Note that the uncertainties in 
the x- and y-directions are used to weight the data in the fit. 



s 2 = 2J (Vi - a ~ bx i) 2 ~ J! ^ + fo2f7 ^) ( 10 ) 

where x and y are abscissa and ordinate values with cor- 
responding uncertainties a x and a y and a and b are the 
y-intercept and slope of of the best fit line. A value of is 
reported if the result from Equation [TU] is negative, as this 
would indicate that all of the scatter around the best fit line 
could be accounted for by the measurement uncertainties. 

Figure [2] shows how the (PAH 8 fim)/2A fim ratio varies 
with 24 ^m surface brightness among the 45 arcsec square 
regions described in Section 12.41 Note that the resolution 
of the data used in this figure is matched to the 38 arcsec 
resolution of the 160 /im images so that the results from 



these figures can be more easily compared to the results in 
Section [4] The slopes and intrinsic scatter for the best fit 
lines are given in Table [3] The best fitting lines are deter- 
mined using uncertainties in both the x- and y-directions to 
weight the data, so the fits are strongly weighted towards 
high surface brightness regions. 

The slopes of the best fit lines in Figure [5] are generally 
not statistically equivalent to 0, which demonstrates that 
the (PAH 8 fim) /24 /im ratio varies with surface brightness. 
In many galaxies (most notably NGC 3031, NGC 3351, and 
NGC 6946), the (PAH 8 /um)/24 /jm ratio decreases notably 
as the 24 /im surface brightness increases. These tend to 
be galaxies with infrared-bright point-like nuclei. In a few 
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Figure 3. Plots of the (PAH 8 /Ltm)/24 fim surface brightness ratios versus deprojected galactocentric radii for the 45 arcsec square 
regions measured in these galaxies. Note that the data are extracted from images with t he same resolution as the 160 fim data. The radii 
are normalized by the radius of the D25 isophote given bv lde Vaucouleurs et all lll99lh . The grey lines are the best fitting lines for the 
relations in each plot; slopes and intrinsic scatter for these fits are given in Table [6] Only the uncertainties in the (PAH 8 fim)/ 160 fim 
ratio are used to weight the data in the fit; the uncertainties in the radius are assumed to be negligible in this analysis. 



other galaxies (NGC 2403 and NGC 7793, for example), the 
(PAH 8 /im)/24 fim ratio increases as the 24 /im surface 
brightness increases. The data for some galaxies also show 
a broad scatter, particularly for galaxies without infrared- 
bright nuclei. For a given 24 fim surface brightness, the 
(PAH 8 /im)/24 fim ratio varies by over a factor of 2 in 
high signal-to-noise regions in NGC 925 and NGC 2403. Al- 
though both of these galaxies are fit with linear relations in 
Figure [2J the broad scatter indicates that such a fit is un- 
realistic, so the (PAH 8) and 24 fim emission must be only 
weakly associated with each other. Also, some data points 
fall well below the best fit lines in Figure [5] particularly in 
the plots for NGC 2403 and NGC 3938. These regions cor- 



respond to very infrared-bright star-forming regions visible 
in Figure [1] 

Figure [1] does not reveal the presence of any obvious 
dependence of the (PAH 8 fim)/24 fim on radius, but it is 
still useful to measure such gradients for comparison with 
abundance gradients, especially since it has been shown 
that the (PAH 8 /xm)/24 fim ratio varies with metallic- 
itv rtEngelbracht et al J 120051: iDale et al.ll2005l; IPraine et all 
120071 ; ICalzetti et al.ll2007l ; lEngelbracht et alj|2008t ). The ~»T 
lation between the (PAH 8 fim)/2A fim ratio and radius 
for the 45 arcsec square regions described in Section 12.41 is 
shown in Figure [3j with slopes for the best fit lines and in- 
trinsic scatter given in Table 2] As in Figure [5] the slopes 
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Table 3. Slopes and Intrinsic Scatter for the Best Fit Line De- 
scribing (PAH 8) / 24 (im Surface Brightness Ratio versus 24 (im 
Surface Brightness Density™ 



Name 


Slope 


Intrinsic 
Scatter 


NGC 


628 


-0.030 ±0.015 


0.20 


NGC 


925 


0.15 ± 0.02 


0.69 


NGC 


2403 


0.090 ± 0.008 


1.91 


NGC 


3031 


-0.270 ±0.011 


0.93 


NGC 


3184 


-0.030 ±0.011 


0.13 


NGC 


3351 


n ^9 -i- n n9 

— U.OZ ZC u.uz 


n 
u 


NGC 


3938 


-0.035 ± 0.015 


0.14 


NGC 


4254 


-0.085 ±0.013 


0.0011 


NGC 


4321 


-0.128 ±0.010 


0.14 


NGC 


4579 


-0.18 ±0.03 





NGC 


4725 


-0.121 ±0.017 


0.12 


NGC 


5055 


-0.041 ±0.007 


0.43 


NGC 


5194 


-0.038 ± 0.008 


0.89 


NGC 


6946 


-0.21 ±0.02 


0.069 


NGC 


7793 


0.043 ± 0.011 


0.13 



a These slopes and intrinsic scatter measurements are for the 
lines fit to the data in Figure [2] The data are measured within 
45 arcsec bins in images with the same resolution as the 160 /im 
images. Intrinsic scatter measurements of are reported when the 
value calculated with Equation \W\ is negative; this indicates that 
measurement uncertainties can account for all of the scatter in 
the relation. 

vary significantly among the galaxies in this sample. Some 
galaxies with infrared-bright nuclei, such as NGC 3351 and 
NGC 4579, have large positive radial gradients in the (PAH 
8 ftm)/ 24 /im ratio, and in NGC 3031 and NGC 6946, the 
infrared-bright nuclei fall below the relation between the 
(PAH 8 /xm)/24 /im ratio and radius. In a few other galaxies 
(most notably NGC 925 and NGC 2403), the radial gradient 
in the (PAH 8 /xm)/24 /xm ratio is negative. As in the re- 
lation between the (PAH 8 /im)/24 /xm ratio versus 24 /jm 
surface brightness, significant scatter is seen in some rela- 
tions between the (PAH 8 /xm) /24 /xm ratio and radius, and 
some infrared-bright regions fall well below the best fit lines, 
as can be seen most clearly in NGC 2403 and NGC 3938. 
Further discussion concerning how these gradients might be 
related to radial gradients in abundances is presented in Sec- 
tion Oy] 

Overall, these data demonstrate that the relation be- 
tween PAH 8 and 24 /xm emission on spatial scales smaller 
than ~ 2 kpc may exhibit a significant amount of scatter. 
Moreover, the relation between the (PAH 8 /xm) /24 /im ra- 
tio and the 24 /xm surface brightness does not vary in a way 
that is easily predictable within these galaxies. 



4 COMPARISONS OF PAH 8 AND 160 /xM EMISSION 

The panels in the right columns of Figure Q] show how 
the (PAH 8 /im)/160 /im surface brightness ratio varies 
within the sample galaxies. The resolution of the maps is 
38 arcsec. As with the (PAH 8 /xm)/24 /xm ratio, the (PAH 
8 /xm)/160 /xm ratio does vary across the discs of these galax- 
ies, although the variations are notably different. The (PAH 
8 /xm)/160 /xm ratio generally appears to increase as the 



Table 4. Slopes and Intrinsic Scatter for the Best Fit Line Describ- 
ing (PAH 8) / 24 /xm Surface Brightness Ratio versus Radius™ 



Name 


Slope 


Intrinsic 
Scatter 


NGC 


628 


-0.07 ±0.02 


0.16 


NGC 


925 


-0.39 ±0.03 


0.23 


NGC 


2403 


-0.34 ±0.02 


1.43 


NGC 


3031 


0.178 ±0.016 


1.64 


NGC 


3184 


0.051 ±0.019 


0.13 


NGC 


3351 


0.59 ±0.05 


0.18 


NGC 


3938 


U.Ul it U.Uo 


U. 10 


NGC 


4254 


0.17 ±0.02 





NGC 


4321 


0.20 ±0.02 


0.17 


NGC 


4579 


0.38 ±0.05 





NGC 


4725 


0.32 ±0.03 


0.040 


NGC 


5055 


0.082 ± 0.017 


0.42 


NGC 


5194 


0.09 ±0.02 


0.85 


NGC 


6946 


0.14 ±0.03 


0.23 


NGC 


7793 


-0.103 ±0.018 


0.091 



a These slopes and intrinsic scatter measurements are for the 
lines fit to the data in Figure [3] The data are measured within 
45 arcsec bins in images with the same resolution as the 160 /xm 
images. Intrinsic scatter measurements of are reported when the 
value calculated with Equation \W\ is negative; this indicates that 
measurement uncertainties can account for all of the scatter in 
the relation. 

160 /xm surface brightness increases. Moreover, the (PAH 
8 /xm)/160 /xm ratio also appears enhanced in large scale 
structures within the discs of these galaxies, such as the spi- 
ral arms in NGC 3031 and NGC 6946. In the centres of a few 
galaxies, the (PAH 8 /xm)/160 /xm ratio appears to decrease 
slightly. This is most apparent in NGC 3184, NGC 4725, 
and NGC 5055. The disc of NGC 4725 contains a ring, so 
the apparent decrease in the central (PAH 8 /xm)/160 /xm 
ratio could partly be related to an enhancement of the (PAH 
8 /xm)/160 /xm ratio in the ring. 

The maps also show that the (PAH 8 /xm)/160 /xm ratio 
does not necessarily peak in the extranuclear star-forming 
regions where local minima in the (PAH 8 /xm)/24 /xm ratio 
are found. This can be seen most dramatically in the im- 
ages of NGC 2403, NGC 3184 and NGC 3938. To show that 
the differences between the (PAH 8 /xm)/24 /xm and (PAH 
8 /xm)/160 /xm ratio maps is not a result of resolution effects 
for these three galaxies, we show both ratio maps for data at 
a resolution of 38 arcsec in Figure[4] Even at this resolution, 
the infrared-bright HII regions located in the northeast part 
of the disc near the nucleus of NGC 2403, in the northeast 
part of the disc in NGC 3184, and in the east part of the disc 
in NGC 3938 are still identifiable in the (PAH 8 /im) /24 /xm 
ratio maps. However, these three infrared-bright HII regions 
are indistinguishable from the nearby dust emission in the 
(PAH 8 /xm)/160 /xm ratio maps. Nonetheless, keep in mind 
that enhancements in the (PAH 8 /xm)/160 /xm ratios are 
still visible in the large scale structures such as the spiral 
arms in NGC 3031 and NGC 6946. 

Figure [5] shows how the (PAH 8 /xm)/160 /xm sur- 
face brightness ratio varies with 160 /xm surface brightness 
among the sample galaxies, and the slopes and intrinsic scat- 
ter for the best fit lines as well as Spearman's correlation co- 
efficients for the data are given in Table [5] Again, the best 
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Figure 4. Images of the (PAH 8 ^tm)/24 (im surface brightness ratio at the resolution of the 160 data (left column), and the (PAH 
8 /im)/160 £tm surface brightness ratio at the resolution of the 160 fira data (right column) for three galaxies with extranuclcar star- 
forming regions that are exceptionally bright at 24 fim. These maps demonstrate how the lower resolution affects the (PAH 8 fim) /24 (im 
ratio maps and how the maps still differ in comparison to the (PAH 8 /im)/160 fim ratio maps. The circle in the top right corner of each 
map shows the 38 arcsec FWHM of the 160 fim beam. See Figure [T] for additional information. 



© RAS, MNRAS 000, [EHS] 



14 G. J. Bendo et al. 



o 



0.02 
0.01 
0.005 

0.02 
0.01 
0.005 

0.02 




| 0.01 

co 

^0.005 



0.02 
0.01 
0.005 



0.005 



NGC 628 ■ 
i i 1 1 1 n | 1 — i i 1 1 1 n | 1 — i i i 1 1 



NGC 3031 -- 



i i i 1 1 n l 1 — i i i 1 1 n l 1 — i i i 1 1 



NGC 3938 -- 
i i i 1 1 n | 1 — i i i 1 1 n | 1 — i i i 1 1 



NGC 4579 -- 
i i i 1 1 n | 1 — i i i 1 1 n | 1 — i i i 1 1 



NGC 5194 -- 



—I I I I I I I I I l_L 



NGC 925 ■ 
i i i i n il 1 — i i i i n il 1 — i i i 1 1 



NGC 3184 -- 



i i i i i n 1 — i i i i i n 1 — i i i 1 1 




NGC 4254 + 
i i i i n il 1 — i i i i n il 1 — i i i 1 1 




NGC 4725 -- 
i i i i n il 1 — i i i i n il 1 — i i i 1 1 




NGC 6946 -- 



I I I I I I l—l— L 



NGC 2403 ■ 
i i i 1 1 n | 1 — i i i 1 1 n | 1 — i i i 1 1 



NGC 3351 



i i i 1 1 n l 1 — i i i 1 1 n l 1 — i i i 1 1 



NGC 4321 ■ 
i i i 1 1 n | 1 — i i i 1 1 n | 1 — i i i 1 1 



NGC 5055 ■ 
i i i 1 1 n | 1 — i i i 1 1 n | 1 — i i i 1 1 



NGC 7793 



_l I I I I I I I L 



10 



100 



10 



100 



10 100 
l v (160 urn) (MJy sr" 1 ) 

Figure 5. Plots of the (PAH 8 £tm)/160 /an surface brightness ratios versus the 160 fim surface brightnesses for the 45 arcsec square 
regions measured in these galaxies. The grey lines are the best fitting lines for the relations in each plot; slopes, intrinsic scatter, and 
correlation coefficients for these fits are given in Table [5] Note that the uncertainties in the x- and y-directions are used to weight the 
data in the fit. 



fitting lines are determined using uncertainties in both the 
x- and y-directions to weight the data. 

For all galaxies in the sample, the (PAH 8 /jm)/160 /im 
ratio generally increases as the 160 /jm surface brightness 
increases, although the slopes of the relations are relatively 
shallow for some galaxies, such as NGC 3031, NGC 3351, 
and NGC 4725. If the slopes of the best fit lines in Fig- 
ure [5] were equivalent to 0, this would indicate that a one- 
to-one correspondence exists between the PAH 8 and 160 /im 
bands. However, since the slopes are instead all positive, this 
indicates that the colours change from low to high surface 
brightness regions. 

The scatter in the data around the best fit lines gen- 
erally appears to be at the 10%-20% level in many cases. 



According to the intrinsic scatter measurement used here, 
the scatter in many of the plots can be explained mostly by 
uncertainties in the measurements. For most galaxies, the 
intrinsic scatter measurements in Table [5] are either simi- 
lar to or notably lower than the values in Table [3] Because 
the data used for Tables [3] and [S] were measured in images 
that were degraded to the resolution of the 160 /im images, 
resolution effects should not be a factor in this comparison. 
Hence, this comparison between the intrinsic scatter mea- 
surements demonstrates quantitatively that the relation be- 
tween PAH 8 and 160 /im emission may exhibit less scatter 
than the relation between PAH 8 and 24 /im emission. 

Also note that very low and very high surface brightness 
45 arcsec regions in NGC 5194 and NGC 5055 fall below the 
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Figure 6. Plots of the (PAH 8 fim)/l60 /an surface brightness ratios versus deprojected galactocentric radii for the 45 arcsec squar e 
regions measured in these galaxies. The radii are normalized by the radius of the D25 isophote given bv lde Vaucouleurs et ail l ll99ll) . 
The grey lines are the best fitting lines for the relations in each plot; slopes and intrinsic scatter for these fits are given in Table [6] Only 
the uncertainties in the (PAH 8 fim)/160 fim ratio are used to weight the data in the fit; the uncertainties in the radius are assumed to 
be negligible in this analysis. 



best fit line. A related phenomenon is visible in NGC 4725, 
where the 45 arcsec regions within the inner ring fall below 
the best fit line in Figure [5] The disparity in the slopes 
between the high and low surface brightness data for some 
galaxies demonstrates that the (PAH 8 /im)/160 /im ratio 
either stops rising or decreases in the high surface brightness 
centres of the galaxies, as can also be seen in the maps of 
the (PAH 8 ^im)/160 /jm ratio in Figure [1] and in the plots 
of the (PAH 8 pim)/160 fim ratio versus radius in Figure [5] 

Figured] illustrates how the (PAH 8 ^m)/160 fim. ratio 
may peak outside the nuclei of nearby galaxies. From the ra- 
tio maps alone, it is apparent that the (PAH 8 pm)/160 fim 
ratio does not necessarily monotonically decrease from the 



nuclei to the edges of the optical discs as was suggested by 
iBendo et ail <|2006h . Both Figured which plots the (PAH 
8 /jm)/160 /jm ratio versus deprojected galactocentric ra- 
dius for 45 arcsec regions in these galaxies, and Table [6J 
which gives the slopes and intrinsic scatter measurements 
for the best fit lines in Figure [6] as well as the Spearman's 
correlation coefficient for the relations, support this conclu- 
sion. First, note that the regions with the highest (PAH 
8 fim)/ 160 fim ratios within some of these galaxies are found 
outside the nucleus. As can be seen in Figure [IJ the regions 
with enhanced (PAH 8 ^m)/160 fim ratios may correspond 
to spiral structure, as is most clearly seen in NGC 3031 and 
NGC 6946. In NGC 4725, the inner ring has the highest 
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Table 5. Results for the Best Fit Line Describing (PAH 
8) / 160 jim Surface Brightness Ratio versus 160 (im Surface 
Brightness" 



Name 


Slope 


Intrinsic 
Scatter 


Spearman's Rank 
Correlation 
Coefficient^ 


NGC 


628 


0.188 ±0.018 





0.85 


NGC 


925 


0.37 ±0.03 





0.79 


NGC 


2403 


0.415 ±0.012 





0.93 


NGC 


3031 


0.11 ±0.02 





0.27 


NGC 


3184 


0.112 ±0.017 


0.034 


0.70 


NGC 


3351 


n ns 4- n no 
U.Uo it u.uz 


U 


u.oy 


NGC 


3938 


0.12 ±0.02 





0.82 


NGC 


4254 


0.093 ± 0.016 


0.013 


0.83 


NGC 


4321 


0.113 ±0.015 





0.78 


NGC 


4579 


0.11 ±0.06 





0.51 


NGC 


4725 


0.07 ±0.03 


0.17 


0.12 


NGC 


5055 


0.111 ±0.011 


0.26 


0.80 


NGC 


5194 


0.159 ±0.012 


0.56 


0.80 


NGC 


6946 


0.179 ±0.013 


0.040 


0.83 


NGC 


7793 


0.300 ±0.019 





0.88 



a These slopes and intrinsic scatter measurements are for the 
lines fit to the data in Figure [5] The data are measured within 
45 arcsec bins in images with the same resolution as the 160 fim 
images. Intrinsic scatter measurements of are reported when the 
value calculated with Equation \W\ is negative; this indicates that 
measurement uncertainties can account for all of the scatter in 
the relation. 

b The Spearman's rank correlation coefficient may have values 
between -1 and 1. A value close to 1 indicates a direct correla- 
tion between two values. A value close to -1 indicates an inverse 
correlation. A value of indicates no correlation. 

(PAH 8 /im)/160 fim ratio, not the nucleus. The intrinsic 
scatter in the fits versus surface brightness is generally lower 
than for those for the fits versus radius. Moreover, the Spear- 
man's correlation coefficients for the relations between the 
(PAH 8 iim)/160 fim ratio and 160 fim surface brightness 
generally has a higher absolute value than the correspond- 
ing correlation coefficients for the relation between the (PAH 
8 /im)/160 fim ratio and radius. These results suggest that 
the ratio may be more strongly dependent on 160 fim surface 
brightness than radius. 

Since the (PAH 8 iim)/160 fim ratio may be dependent 
on dust heating, we also examine how the ratio is related 
to the 24 fim/160 fim ratio. The 24 /im band, which traces 
>100 K hot dust emission, increases faster than other in- 
frared band (including the PAH 8 and 1 60 fim bands) as 
the illuminating r a diation field increa ses IjDale et alj|200lt 
iLi fc Drainel 1200 ll ; iDraine fc Li I2007T ). The 160 /im band 
is approximately directly proportional to the TIR flux, as 
discussed in Section [2.11 Therefore, the 24 /im/160 /jm ra- 
tio should be a reasonable indicator of dust heating. If the 
mass fraction of PAHs remains constant and if the (PAH 
8 /im)/160 fim ratio is dependent on dust heating, then the 
(PAH 8 fim)/ 160 fim ratio should monotonically increase 
as the 24 /im/160 fim ratio increases, although the slope 
may not necessarily be constant. This comparison is simi- 
lar to the direct comparison between the PAH 8 and 24 fim 
bands performed in Section [3] but the normalisation with 
the 160 fim band removes variations related to dust surface 
density. 



Table 6. Slopes and Intrinsic Scatter for the Best Fit Line Describ- 
ing (PAH 8) / 160 fim Surface Brightness Ratio versus Radius™ 



Name 


Slope 


Intrinsic 
Scatter 


Spearman's Rank 
Correlation 
Coefficient 6 


NGC 628 


-0.21 ±0.03 


0.036 


-0.80 


NGC 925 


-0.42 ±0.04 


0.038 


-0.74 


NGC 2403 


-0.83 ±0.03 


0.29 


-0.88 


NGC 3031 


0.05 ±0.02 





0.11 


NGC 3184 


-0.125 ±0.03 


0.11 


-0.68 


NGC 3351 


—0.12 ± 0.04 





-0.68 


NGC 3938 


-0.18 ±0.04 





-0.84 


NGC 4254 


-0.11 ±0.02 


0.039 


-0.74 


NGC 4321 


-0.16 ±0.02 


0.00177 


-0.73 


NGC 4579 


-0.04 ±0.08 





-0.32 


NGC 4725 


0.12 ±0.05 


0.130 


0.18 


NGC 5055 


-0.23 ±0.02 


0.28 


-0.79 


NGC 5194 


-0.22 ±0.02 


1.07 


0.71 


NGC 6946 


-0.188 ±0.016 


0.069 


-0.73 


NGC 7793 


-0.42 ±0.03 





-0.90 



a These slopes and intrinsic scatter measurements are for the 
lines fit to the data in Figure [6] The data are measured within 
45 arcsec bins in images with the same resolution as the 160 fim 
images. Intrinsic scatter measurements of are reported when the 
value calculated with Equation \W\ is negative; this indicates that 
measurement uncertainties can account for all of the scatter in 
the relation. 

b The Spearman's rank correlation coefficient may have values 
between -1 and 1. A value close to 1 indicates a direct correla- 
tion between two values. A value close to -1 indicates an inverse 
correlation. A value of indicates no correlation. 

Figure [7J shows how the (PAH 8 /xm)/160 /im ratio 
varies with the 24 /xm/160 /im ratio within 45 arcsec re- 
gions in the sample galaxies. Many of the 45 arcsec regions 
that are relatively weak in 24 fim emission tend to show 
a tight correspondence between the (PAH 8 fim)/160 fim 
and 24 /im/160 fim ratios, but some regions with enhanced 
24 fim emission appear far to the right of these curves. The 
most spectacular example is NGC 3938, where the ratios for 
most of the regions closely follow a linear relation but the 
HII region on the east side of the disc falls far to the right 
of this relation. The relation between the two ratios on the 
left side of the panels in Figure [7] suggest that the (PAH 
8 /im)/160 fim ratio is dependent on dust heating, but the 
outliers on the right in these panels show PAH 8 fim emis- 
sion is not enhanced in areas with strong dust heating such 
as HII regions. Further interpretation of these results is pre- 
sented in the next section. 

For another perspective on the nature of the variation 
in the (PAH 8 /xm)/160 fim ratio, we examined how the 
160 /xm/TIR and (PAH 8 fim) /TIR ratios vary as a function 
of TIR surface brightness. These are displayed in Figures [8] 
and [9] with slopes for the best fit lines given in Table [7J 
and Table As stated in Section \TT\ the 160 /xm/TIR ra- 
tio should gradually decrease as dust temperatures increase 
above ~ 15 K. This is reflected in the general trends vis- 
ible in Figure [5] Some scatter is visible in this figure, but 
this is expected, as the ratio will deviate from this general 
trend where 24 fim emission is strongly enhanced within star 
forming regions. The (PAH 8 /xm)/TIR ratio is constant or 
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Figure 7. Plots of the (PAH 8 /xm)/160 jim surface brightness ratios versus the 24 /im/160 (im surface brightness ratios for the 45 arcsec 
square regions measured in these galaxies. The 24 fim/160 fim ratios are used as a proxy for dust heating here. 



changes very little in many galaxies (e.g. the slope of the best 
fit line in Figure deviates less than 3<r from 0), which indi- 
cates that the (PAH 8/im) emission is directly proportional 
to TIR emission in these galaxies. Some scatter related to the 
scatter seen in Figure [5] is visible in these data, and in some 
cases, the infrared-brightest regions near the nucleus deviate 
from the trend in (PAH 8 /im) /TIR ratio versus TIR surface 
brightness, which makes the best fit lines appear positive. In 
a few exceptional galaxies, however, the (PAH 8 pm)/TIR 
ratio increases as TIR surface brightness increases. We con- 
clude that the variations in (PAH 8 ^m)/160 /im ratio are 
driven by changes in 160 /jm emission relative to total dust 



emission in all galaxies and also by changes in (PAH 8 fim) 
emission relative to total dust emission in some situations. 
This is discussed further in the following section. 



5 DISCUSSION 

5.1 Interpretation of the relation between PAH 8 and 
24 fj,m emission 

The results in Section [3] as well as [4] demonstrate that the 
relation between PAH 8 and 24 fim emission exhibits a sig- 
nificant scatter on spatial scales smaller than ~ 2 kpc. In 
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Figure 7 - continued 



contrast, some ISO results from comparisons of 7 and 15 /im 
emission had implied that the ratio of PAH to hot dust emis- 
sion should be relatively uniform across the discs of most 
galaxies (e.g. iRoussel et al]|200ll ). However, some ISO stud- 
ies ac tually found sign ificant variations in the 7 /im/15 /im 
ratio (Haas et al. 2 0021 ') . and even studies that did find uni- 
form 7 /xm/15 /jm colours within the discs of galaxies noted 
that the ratio might d ecrease within starb ursts in the cen- 
tres of some galaxies (|Roussel et al.ll200ll ). The variations 
in the (PAH 8 /im)/24 /im ratio observed in this sub- 
sample of SINGS galaxies are consistent wi th Spitzer re- 
sults for indi vidual galaxies, such as for M 51 l|Calzetti et al.l 
120051), M81 l|Perez-Gonzalez et all 12001 ) . and NGC 4631 
l|Bendo et al.ll2006T) . Moreover, we find that 24 /im emis- 
sion is more peaked in the centres of HII regions while the 
PAH 8 /mi emission is relatively stronger outside HII re- 
gions, which i s consistent with s imilar phenom ena observed 
in N GC 300 (iHelou et al.ll2004[), NG C 4631 l|Bendo et all 
120061 ), and M101 (jGordon et alj|2008h . 

The variations in the (PAH 8 /mi) /24 /im ratio are best 
seen by comparing bright star-forming regions and diffuse 
regions. If, because of the limited sensitivity of the data, 
the ratio between PAH and hot dust emission is only mea- 
sur ed in the bright reg ions in some galaxies, as was done 
by IRoussel et al] (1200 ll ). then the ratio between PAH and 
hot dust emission may appear uniform. When these bright 
regions are compared to diffuse emission, however, varia- 
tions may be seen in the ratio of PAH to hot dust emission. 



This may be the primary reason why some results from ISO 
suggested that the ratio of PAH to hot dust emission was 
uniform in most spiral galaxies whereas variations may be 
seen in Spitzer data. 

Additionally, differences between the nature of 15 and 
24 fim dust emission could have also contributed to the dif- 
ferences between the 7 /im/15 /im relation observed with 
ISO and the (PAH 8 /im)/24 /im relation observed with 
Spitzer. Based on semi-empirical and theoretical models, the 
24 /im band is expected to increase faster t han the 15 /im as 
the illuminating radiation field incre ases (|Dale et al.ll200l1 ; 
iLi fc Drainell200ll ; iDraine fc Lill2007l ). Moreover, unlike the 
24 /im band, a signif icant fraction of the 15 /im band in - 
cludes PAH emission l|Smith et al.ll2007l ; [Praine fc Lill2007f l. 
Consequently, a comparison of PAH to 24 /im data should 
exhibit more scatter than a comparison of PAH to 15 /im 
data. 

We also found that differences in the distribution of 
star-forming regions within galaxies could influence how the 
(PAH 8 /im)/24 /im ratio varies with 24 /im surface bright- 
ness. In some galaxies, such as NGC 3351 and NGC 6946, 
the nucleus is the strongest site of star formation. Hence, the 
location expected to have the lowest (PAH 8 /im) /24 /im ra- 
tio will correspond to the location with the highest 24 /im 
surface brightness, so the data will show that the (PAH 
8 /im)/24 /im ratio decreases as the 24 /im surface bright- 
ness increases, thus va rying as expected from many models 
of dust emission (e.g. iDale et alJl200ll ; ILi fc Drainej 1200 ll ; 
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Figure 8. Plots of the 160 /im/TIR surface brightness ratios versus TIR surface brightness for the 45 arcsec square regions measured in 
these galaxies. The grey lines are the best fitting lines for the relations in each plot; slopes for these fits are given in Table [7] Note that 
the uncertainties in the x- and y-directions are used to weight the data in the fit. 



iDraine fc Lill2007l ). In other galaxies, such as NGC 925 and 
NGC 2403, the nucleus is not the site of the strongest star 
formation activity, and many HII regions can be found at the 
periphery of the regions that were detected at the 3a level 
in these data. While p oint-like 24 /xm sources should cor- 
respo nd to HII regions (|Prescott et al ] |2007l ; ICalzettT et al.l 
20071). dust e missio n models su c h as those presented by 



Li fc Drainel (|200lh . bale et al.1 (120011 ). and IDraine fc l! 
1 20071 ) suggest that diffuse 24 /im dust emission may still 



potentially originate from regions outside star-forming re- 
gions with high radiation fields. Furthermore, ionising pho- 
tons may escape from s tar-forming reg ions into the dif- 
fuse interstellar med ium (|Oev et al.ll2007| ). and results from 
ICalzetti et al.l |2007l ) suggest that 24 fj,m emission may trace 
the recombination line emission from the diffuse gas heated 



by these photons. Hence, it is possible that the infrared- 
brightest regions in the centres of some galaxies may contain 
a higher fraction of 24 /im emission from diffuse dust than 
the infrared-faint regions, so the (PAH 8 /im)/24 fim ratio 
will increase as the 24 surface brightness increases, as 
can be seen for NGC 925 and NGC 2403 in Figured 

Several mechanisms may be responsible for creating 
the differences in the spatial distribution of the PAH 8 
and 24 /im emission. First, as the illuminating radiation 
field increases, the 24 /im band is thought to inc rease more 
rapidly than bands that trace PAH emis sion l|Dale et al.l 
l200ll ; (Li fc Drainel l200ll ; IDraine fc Lill2007l ). Therefore, the 
24 fim emission may be more strongly enhanced than the 
PAH emission in regions with very high radiation fields, such 
as star-forming regions. Second, the PAHs may be destroyed 
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Figure 9. Plots of the (PAH 8 fim) /TIR surface brightness ratios versus TIR surface brightness for the 45 arcsec square regions measured 
in these galaxies. The grey lines are the best fitting lines for the relations in each plot; slopes for these fits are given in Table [8] Note 
that the uncertainties in the x- and y-directions are used to weight the data in the fit. 



in regions with strong radiation fields o r in regions with larg e 
numbers of high-energy photons (e.g. iMadden et "aTl l2006). 
The PAHs would be absent, and the (PAH 8 /im)/24 /xm 
ratio would be low in the centres of star-forming regions 
and AGN. Third, the ratio of the 7.7 /an PAH feature to 
other PAH features may vary with changes in the charge 
state of the PAHs (e.g. | Allamandola et~al .1 19991 : iLi fc Draind 
l200ll ; iDraine fc Liil2007h . This could occur in HII regions if 
the electron densities are high enough that recombination 
rates becomes significant compared to photoionisation rates, 
which would lead to a redu ction of the PAH + ions that pro- 
duce the 7.7 /jm emission (|Weingartner fc Draind liOQlh . 

While the enhancement of 24 /im emission linked to dust 
heating in star-forming regions is certainly at least part of 
the reason why the (PAH 8 /im)/24 /mi ratio varies, we ar- 



gue that the variations must also be caused in part by a 
reduction in PAH emission in the 8 /im band through ei- 
ther PAH destruction or changes in the relative strengths 
of PAH spectral features. Models of dust emission have sug- 
gested that, if PAHs are present and radiating at 8 /zm, 
the (PAH 8 /im)/160 /im ratio would be enhanced in re- 
gions with very strong radiat ion fields, warmer dust, and 
low (PAH 8 /im) /24 /im ratios (|Draine fc Lill2007T I. However, 
the (PAH 8 /im)/160 /im ratio does not peak within indi- 
vidual star-forming regions where local minima in the (PAH 
8 /im)/24 fim ratio are found, such as the bright extranu- 
clear regions in NGC 2403, NGC 3184, and NGC 3938. This 
is best demonstrated by comparisons of the maps of (PAH 
8 /im)/24 /zm and (PAH 8 /mi)/160 /mi ratios at matching 
resolutions in Figure [3] although the additional discussion 
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Table 7. Results for the Best Fit Line Describing 160 /mi / TIR 
Surface Brightness Ratio versus TIR Surface Brightness 



Name 



Slope" 



NGC 


628 


—0.092 


± 





011 


NGC 


925 


—0.098 


± 





017 


NGC 


2403 


—0.130 


± 





005 


NGC 


3031 


—0.186 


± 





012 


NGC 


3184 


-0.057 


± 





009 


NGC 


3351 


-0.173 


± 





013 


NGC 


3938 


-0.070 


± 





015 


NGC 


4254 


-0.088 


± 





014 


NGC 


4321 


-0.113 


± 





009 


NGC 


4579 


-0.15 


± 





02 


NGC 


4725 


-0.083 


± 





013 


NGC 


5055 


-0.089 


± 





006 


NGC 


5194 


-0.082 


± 





008 


NGC 


6946 


-0.138 


± 





013 


NGC 


7793 


-0.094 


± 





009 



a These slopes are for the lines fit to the data in Figure [8] The 
data are measured within 45 arcsec bins in images with the same 
resolution as the 160 /urn images. 



Table 8. Results for the Best Fit Line Describing (PAH 
8 (im) / TIR Surface Brightness Ratio versus TIR Surface Bright- 
ness 



Name 


Slope" 


NGC 


628 


0.073 ±0.012 


NGC 


925 


0.237 ±0.018 


NGC 


2403 


0.210 ±0.006 


NGC 


3031 


-0.086 ±0.011 


NGC 


3184 


0.037 ± 0.009 


NGC 


3351 


-0.087 ±0.012 


NGC 


3938 


0.048 ±0.012 


NGC 


4254 


0.000 ± 0.010 


NGC 


4321 


-0.018 ± 0.007 


NGC 


4579 


-0.06 ± 0.03 


NGC 


4725 


-0.028 ±0.013 


NGC 


5055 


0.013 ±0.006 


NGC 


5194 


0.062 ± 0.006 


NGC 


6946 


0.024 ± 0.008 


NGC 


7793 


0.183 ±0.010 



a These slopes are for the lines fit to the data in Figure [9] The 
data are measured within 45 arcsec bins in images with the same 
resolution as the 160 /tm images. 



in Section [3] and the maps in Figure [T] also support this con- 
clusion. Furthermore, Figure [7] demonstrates that the (PAH 
8 /mi)/ 160 /jm ratios are not always higher in regions with 
stron gly enhanced 24 [im emission. According to the mod- 
els of iDraine fell (|2007m . the (PAH 8 /mT)/160 /im should 
monotonically increase with the 24 /mi/ 160 /mi ratio in these 
regions if PAHs are present and if the PAH ionisation does 
not change. These results indicate that the variations in the 
(PAH 8 /mi)/24 /mi ratio must be in part caused by the 
suppression of PAH 8 /mi emission. 

Observations of individual HII regions within the 
Milky Way have shown that the PAH emission may be 
found primarily in shell-like structures around the star- 



forming regions (e.g. Rho et al 1l2006l ; IChurchwell et al ]|2006l ; 
ISmith fc Brooksi l2007h and that the strength of all PAH 
spectral features relative to 24 /mi hot dust emission de- 
creases within the centres of HII regions (e.g. iPovich et al.l 
120071 : iLebouteiller et al.i r2007). If these results are applica- 
ble to the galaxies in this paper, then the inferred decrease 
in the (PAH 8 /im) /24 /mi ratio within star-forming regions 
may also be partly caused by PAH destruction and not just 
variations in the strengths of PAH spectral features. Further 
work with mid-infrared spectroscopic observations of star- 
forming regions within these galaxies, specifically studies of 
how the total emission in all PAH spectral feature emission 
varies with respect to total dust emission, would be needed 
to confirm that PAH destruction is taking place. 

Although the PAH 8 /im emission does not share a one- 
to-one correspondence with 24 /im emission or with other 
star formation tracers on kiloparsec scales, this does not 
necessarily preclude its use as a tracer of integrated star for- 
mation within nearby spiral galaxies. The variations between 
PAH 8 /im emission and other star formation tracers could 
potentially be averaged out when integrating across the op- 
tical discs of spiral galaxies. However, because a significant 
fraction of the PAH 8 /im emission still originates from the 
diffuse ISM, it would not be as reliable as other global star 
formation tracers. Furthermore, because the ratio of PAH to 
TIR flux and the ratio of PAH to 24 flux density varies 
with metallicity jEngelbracht et~al] 120051: iDale et alj|2005l: 
Draine et all 120071 ; ICalzetti et all 120071 ; lEngelbracht et al l 



200l), a conv^ion factor relating— egrated PAH 8 ,m 



emission to star formation may also vary with global metal- 
licity. 

We also note that variations in metallicity could pos- 
sibly produce variations in the (PAH 8 /mi)/24 /im ratio 
observed in the galaxies studied in this paper, but these 
variations would not be expected to cause the differences in 
the ratio between star-forming regions and the diffuse ISM 
observed here. Moreover, the gas-phase abundances of the 
galaxies in our sample are well above 12±log(0/H)~8-8.1, 
where P AH 8 emi s sion s houl d be suppressed as deter - 
mined bv lDraine et al.l (|2007T ) and lEngelbracht et all l|2008h . 
although NGC 628, NGC 925, NGC 2403, and NGC 3031 all 
have characteristic abundances close to this value and might 
be affected. To examine this issue further, we compared the 
radial gradients in the (PAH 8 /mi)/24 /mi ratio to the radial 
gradients in 12±log(0/H) measu red by J. Moustaka s et al . 
(2008, in preparation) using the IPilvugin fc Thuanl (|2005h 
calibration for 13 of the galaxies in our sample. (Abundance 
gradients are not given by J. Moustakas et al. for NGC 3938 
and NGC 4579, and so these galaxies are not included in this 
comparison.) The relation between these two radial gradi- 
ents can be seen in Figure 1101 If the variations in the (PAH 
8 /mi)/24 /im ratio depended primarily on metallicity vari- 
ations, then we would expect to see a correlation between 
the gradients in Figure 1101 However, the two gradients ap- 
pear uncorrelated, which indicates that the variations in the 
(PAH 8 /im) /24 /im ratio that have been observed here de- 
pend primarily on effects unrelated to metallicity. This con- 
clusion for these spiral galaxies is also consiste nt with the re- 
sults p resented for the spiral galaxy M101 by I Gordon et al.l 
(2008), who found that PAH equivalent widths appeared 
more dependent on the ionisation of the ISM (characterised 
using [NeIII]/[NeII] and [SIV]/[SIII]) than 12+log(0/H). 
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Figure 10. Plot of the radial gradients in the (PAH 8 /tm)/24 (im 
ratio from Table [4] versus radial gradients in 12+log(0/H) mea- 
sured by J. Mous t akas et al. (2008, in preparation) using the 
iPilvugin fe Thuanl (2005) calibration. The gradients are in units 
of dex divided by the radius of the D25 isophote. Gradients in 
12+log(0/H) were not provided for NGC 3938 and NGC 4579. 



5.2 Interpretation of the relation between PAH 8 and 
160 /im emission 

The (PAH 8 /mi)/160 /im ratio appears to be closely cor- 
related with the 160 /mi surface brightness measured on 
45 arcsec scales in most galaxies, even in nearby galaxies 
where this angular scale corresponds to <1 kpc. Moreover, 
the (PAH 8 /tm)/160 /mi ratio sometimes traces high surface 
brightness large scale structures in the discs of these galax- 
ies, such as the spiral arms in NGC 3031 or NGC 6946. This 
indicates that the variations in the ratio may not be primar - 
ily dependent on radius as inferred by iBendo et al.l (|2006t ) 
but instead may be primarily dependent on the 160 /mi sur- 
face brightness. However, the large scale structures within 
these galaxies are only marginally resolved in the 160 fim 
data. Higher resolution observations at wavelengths longer 
than 100 /im are needed to confirm that this interpretation 
is valid. 

Nonetheless, if the above interpretation is correct, 
then it also suggests that the variations in the (PAH 
8 /zm)/160 /jm ratio may be more dependent on 160 /mi 
surface brightness than metallicity within the regions 
studied in these galaxies. Since variations in metallicity 
have been linked to decreased PAH emission relative to 
longer-wavelengt h dust emission in the i ntegrated spectra 
of galaxies (e.g. lEngelbracht et all 120051; [Dale et al.1 1200I 



iDraine et ail 12007^ lEngelbracht et al.l I2OO8) . it was not un- 
reasonable to expect that the observed variations in the 
(PAH 8 /mi)/ 160 /im ratio within these galaxies might be 
linked to metallicity variations. However, the 12+log(0/H) 
values measured in these galaxies by J. Moustakas et al. 
(2008, in prepa ration) genera l ly do not d rop below ~ 8-8.1, 
which is where IDraine et al.l |2007l ) and lEngelbracht et al] 
(2008) showed that metallicity strongly affects PAH 8 /im 



Figure 11. Plot of the radial gradients in the (PAH 8 /tm)/160 /tin 
ratio from Table [6] versus radial gradients in 12+log(0/H) mea- 
sured by J. Mous t akas et al. (2008, in preparation) using the 
IPilvugin fc Thuanl <2005l) calibration. The gradients are in units 
of dex divided by the radius of the D25 isophote. Gradients in 
12+log(0/H) were not provided for NGC 3938 and NGC 4579, 
and the gradients in the (PAH 8 /mi)/160 (im ratio for NGC 4725 
changes significantly between the centre and outer disc, so these 
three galaxies are not included in this plot. 



emission. Moreover, metallicity is expected to decrease 
monotonically with radius, while the (PAH 8 /mi)/160 fim 
ratio and the 160 /im surface brightness do not, and the 
metallicity should not peak within substructures such as 
spiral arms, whereas the (PAH 8 /tm)/160 /mi ratio and 
the 160 /tm surface brightness both peak within such sub- 
structures. As an additional test, we compared the gradi- 
ents in the (PAH 8 /tm)/160 /im ratio versus radius with 
the metallicity gradients from J. Moustakas et al. (2008, in 
preparation) calculated with the IPilvugin fc Thuanl (|2005h 
calibration. For this comparison, we excluded NGC 4725 be- 
cause the gradients in the (PAH 8 /tm)/160 /im ratio changes 
significantly between the nucleus and the outer disc, and we 
excluded NGC 3938 and NGC 4579 because abundance gra- 
dients are not given by J. Moustakas et al. The gradients 
for the other galaxies are plotted in Figure 1111 If the (PAH 
8 /tm)/160 /im ratio was affected by metallicity in the re- 
gions studied in these galaxies, then these data would be 
positively correlated. Since the data in Figure [TT] do not 
exhibit such a correlation, the two gradients may be unre- 
lated. We therefore conclude that, in the regions of the galax- 
ies studied here, metallicity variations are not as important 
as 160 (im surface brightness variations in determining the 
(PAH 8 /tm)/160 /tm ratio, although metallicity may be a 
factor outside the optical d iscs. Again, thi s is con sistent with 
the conclusions reached by I Gordon et al.l l|2008l ). who found 
that PAH equivalent widths in M101 were more dependent 
on the ionisation of the ISM than abundances. 

While the (PAH 8 /mi)/160 /im ratio traces large scale 
structure, we have demonstrated that the ratio is not en- 
hanced within individual star-forming regions, and we ex- 
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plained earlier in this section that PAH 8 fim emission must 
be inhibited in regions with strong radiation fields. Based 
on these conclusions and the strong relation between the 
(PAH 8 fim)/ 160 fim ratio and 160 fim surface brightness, 
we infer that the PAHs in these galaxies are generally as- 
sociated with the cold (~ 20 K) dust that dominates the 
160 fim emission, at least on scales of ~ 2 kpc. Because 
most of this cold dust may be expected to be found in the 
diffuse ISM, the PAHs may also be found primarily in the 
diffuse ISM as well, although some of the cold dust and 
PAHs may also be found within clouds associated with star- 
forming regions. Moreover, since the (PAH 8 ftm)/160 fim 
ratio increases as the 160 fim surface brightness increases, 
the (PAH 8 ftm)/160 /jm ratio may be an indicator of varia- 
tions in the intensity of the radiation field heating the diffuse 
ISM. This is further supported by the tight correlation be- 
tween the (PAH 8 fim)/160 fim ratio and the 24 fim/160 fim 
ratio found for many regions with weak 24 fim emission in 
Figure [7] which are presumably regions that primarily sam- 
ple dust emission from the diffuse ISM. Nonetheless, this 
interpretation is only valid if the PAH mass fraction does 
not vary appreciably between infrared-faint and infrared- 
bright regions in the diffuse ISM, and far-infrared observa- 
tions with higher angular resolution will be needed to deter- 
mine whether this association is still applicable on smaller 
spatial scales. 

In most galaxies in this sample, the variation in the 
(PAH 8 /im)/160 fim ratio with 160 fim surface brightness 
is purely driven by the decrease in 160 /im emission rel- 
ative to TIR emission, as the (PAH 8 fim) /TIR ratio re- 
mains constant with TIR surface brightness. In a few excep- 
tions, however, a decline in the PAH 8 /jm emission relative 
to TIR emission is also partly responsible for the observed 
variations. This reduction of PAH emission may occur if 
the fraction of starlight from evolved red stars increases 
as the 160 /im surface brightness decreases. In this sce- 
nario, the photons in the illuminating radiation field would 
have reduced energies, and the peak temperatures attained 
by PAHs following the absorpiton of single photons would 
be reduced. Consequently, the PAHs would tend to radiate 
at longer wavelengths, and the (PAH 8 /im)/160 /im ratio 
would appear to decrease. 

We also noted that the (PAH 8 fim)/160 fim ratio drops 
within the high surface brightness centre of some galaxies 
within this sample. A couple of mechanisms could be re- 
sponsible for this phenomenon. First, AGN may be respon- 
sible for inhibiting the PAH 8 fim emission either by chang- 
ing the ionisation state of the P AHs or by de s troyin g the 
PAHs. This was also suggested bv lSmith et al.1 (|2007T l. who 
showed that the 7.7 fim PAH spectral feature was suppressed 
within the centres of SINGS galaxies with low luminosity 
AGN. Many of the galaxies we observed with suppressed 
(PAH 8 fim)/160 fim ratios in their centres are objects that 
are classified as containing AGN, although NGC 3184 is 
a notable exception. Another possibility is that the (PAH 
8 fim)/ 160 fim ratio drops in regions where a significant frac- 
tion of the interstellar radiation field originates from evolved 
bulge stars. As explained above, PAHs may not be heated 
as strongly by single photons in such radiation fields, and 
therefore 7.7 fim emission will be inhibited. 



6 CONCLUSION 

We have shown that the relation between PAH 8 and 24 fim 
emission in this sample of spiral galaxies exhibit significant 
scatter on scales of ~ 2 kpc. In particular, we have shown 
that the PAH 8 fim emission is relatively weak compared 
to 24 fim emission in star-forming regions, but the PAH 
8 fim emission is relatively strong in the diffuse ISM. In some 
cases, the (PAH 8 /im)/24 fim ratio may exhibit variations 
greater than a factor of 2 at a given 24 fim surface brightness. 
We argue that two mechanisms are responsible for these 
variations: enhancements in the 24 fim emission relative to 
PAH emission in regions with very strong radiation fields 
with high-energy photons and the corresponding reduction 
of PAH 8 fim emission. 

We have also shown that the PAH 8 fim emission is as- 
sociated with 160 fim emission on scales of ~ 2 kpc within 
these spiral galaxies and that the (PAH 8 fim)/ 160 fim ratio 
appears to be a function of 160 fim surface brightness. The 
scatter in the relation is only at the 10%-20% level for a 
given surface brightness, and the intrinsic scatter measured 
in the relation between the (PAH 8 fim)/ 160 fim ratio and 
160 fim surface brightness is notably lower than for the re- 
lation between the (PAH 8 fim)/24 fim ratio and 24 fim 
surface brightness. While the (PAH 8 fim) / 160 fim ratio ap- 
pears to decrease monotonically with radius in some galax- 
ies, the presence of peaks in the (PAH 8 iim)/160 fim ratio 
corresponding to large scale structure indicates that surface 
brightness may be a more important factor, and statistical 
tests of the best fit lines between the (PAH 8 fim)/ 160 fim 
ratio and either 160 fim surface brightness or radius showed 
that the dependence with 160 fim surface brightness was 
more significant. The strong correlation between the (PAH 
8 fim)/ 160 fim ratio and 160 fim emission and the results 
from the comparison of PAH 8 and 24 fim emission suggest 
that most of the PAHs are located in the diffuse interstellar 
medium with the dust grains that produce the majority of 
the 160 fim emission. We therefore suggest that the (PAH 
8 fim)/ 160 fim ratio may be indicative of the intensity of the 
interstellar radiation field that heats the diffuse interstellar 
dust in these galaxies. 

The results here indicate that PAH emission should be 
used very cautiously as a tracer of star formation on kilo- 
parsec scales if it should be used at all. Instead, the PAH 
emission may be more indicative of the distribution of dif- 
fuse dust within nearby galaxies, although PAH emission 
still may b e affect e d by metallicity (lEngelbracht et alj|200a: 
Dale et all 120051 ; ICalzetti et~aiT 120071 ; lEngelbracht et al I 
2008) , and the results here suggest that PAH emission might 
be inhibited in the diffuse ISM if the radiation field is very 
high. Follow-up observations with the Herschel Space Ob- 
servatory and the James Clerk Maxwell Telescope will allow 
for studying the correlation between PAH and ~ 20 K dust 
emission on smaller spatial scales in nearby galaxies, thus 
placing further constraints on the relation between PAHs 
and cool dust. 
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